Magnetic Seed Field Generation from Electroweak Bubble Collisions, with Bubble 

Walls of Finite Thickness 
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Building on earlier work, we develop an equation-of-motion method for calculating magnetic seed 
fields generated from currents arising from charged W fields in bubble collisions during a first- 
order primordial electroweak phase transition allowed in some proposed extensions of the Standard 
Model. The novel feature of our work is that it takes into account, for the first time, the dynamics of 
the bubble walls in such collisions. We conclude that for bubbles with sufficiently thin surfaces the 
magnetic seed fields may be comparable to, or larger than, those found in earlier work. Thus, our 
results strengthen the conclusions of previous studies that cosmic magnetic fields observed today 
may originate from seeds created during the electroweak phase transition, and consequently that 
these fields may offer a clue relevant to extensions of the Standard Model. 
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PACS Indices:12.38.Lg,12.38.Mh,98.80.Cq,98.80Hw 



INTRODUCTION 

Explaining the origin of galactic and extra-galactic 
magnetic fields remains an outstanding problem in cos- 
mology. A common approach is to view these fields as 
having arisen from magnetic seed fields created in the 
early universe, possibly during the electroweak phase 
transition. Several papers 0, 0,13] have shown how mag- 
netic fields can arise from the equilibration of the Higgs 
phase within the collision region of the expanding elec- 
troweak bubbles in first-order phase transitions. First- 
order phase transitions are notpossible for the Standard 
Model Q but may be allowed [!, H, 0] in certain minimal 
extensions of the Standard Model (MSSM). 

These early studies of magnetic seed formulation were 
formulated in the context of a simple Abelian model. In 
Q we have proposed an alternative and, we believe, a 
more natural mechanism, which is that the charged W- 
fields are the physical origin of the electromagnetic cur- 
rents creating the magnetic seed fields. We found that 
the magnetic seed fields generated by this mechanism 
were of a magnitude comparable to those found in the 
Abelian model. 

In Q we obtain our results by solving equations of 
motion (EOM) for the charged IT-fields of the MSSM 
for bubbles that have collided, 
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d v d ■ w a + m 2 w° = , 



(1) 



where a — (1,2) are the fields of the charged W gauge 
fields and m 2 = p\g 2 jl is the square of its mass within 
the bubbles. Our method of solution followed previous 
studies with jump boundary conditions imposed on the 
w z field at the moment of collision. In this approach, 
the divergence d ■ w a vanishes, a result that follows from 



taking the divergence of both sides of (JTJ) , 
d v m 2 w v = m 2 d u w v = 



(2) 



and taking m(x) to be constant within the bubbles. 

In the present work, we employ a generalization [9l[ic| 
of Ref. [8| that evolves the collision from bubbles that 
are initially separated to estimate, using representative 
boundary conditions, the importance of bubble surface 
dynamics on magnetic field creation. In this case, m 
can no longer be taken constant in x, and the divergence 
of w is now expressed in terms of m(x) as an auxiliary 
condition, 



d-w + X;w v d u m 2 = , 
m 



(3) 



which the solutions of the EOM must the divergence of 
w u in the EOM represents a coupling to the bubble wall 
that cannot be neglected as a contribution to the mag- 
netic field. When surface is taken into account, solving 
Eq. (JTJ subject to Eq. (J3J) requires a new approach, and 
the main difficulty in implementing it is having to solve 
the field equations numerically, in contrast to Ref. 
We note in passing that for the case of bubbles with in- 
finitely thin walls, dm/dr becomes a delta function at 
the bubble surface, and consequently the VF-fields (and 
hence the magnetic fields) may become singular in this 
limit 
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Once the EOM have been solved for the W fields, the 
electromagnetic current may be calculated [1, 0, [l(| from 
the expression 

4nj v = Ge ab3 (w b v {d ■ w a ) - w%d v w» h + 2w a ^w b v ) (4) 

where 
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As in d, 0, the current will vanish if w h v (x) oc 
To ensure a non- vanishing current in the present 
work, we choose boundary conditions at a time t = t = 
such that the fields w z (t ,x) = w z (t ,x) in one bub- 
ble, referred to as boundary condition II (BCII), and 
w z (to,x) = —w®(to,x) in the other (BCI) with the val- 
ues of w u (t ,x) constant in each bubble as in Ref. Q. 
Notice that fixing the boundary condition on w z in this 
way follows Ref. [8(, where the u^-field was taken to be a 
step function at to- Initial conditions for the other fields 
needed in our present work that result from these bound- 
ary conditions and the auxiliary condition are discussed 
below. 

Other choices of boundary conditions for w z may also 
be envisioned, and in principle all these should be aver- 
aged over to find the magnetic field. However, BCI and 
BCII are representative boundary consitions [8| and are 
thus sufficient for our present purposes. 



BUBBLE COLLISION IN CYLINDRICAL 
COORDINATES 

To study the case of two colliding bubbles, we will use 
the axial symmetry about the x\ and X2 axes to write 
the w vector fields in cylindrical coordinates as: 
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In cylindrical coordinates, the auxiliary condition be- 
comes: 



d ■ w(x) 
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Then the equations for the W-fields may then be written 



d 2 wo d 2 wo 



dt 2 



d 2 w 



1 dwo 
r dr 



dr 2 
d 

—d ■ w(x) = 
dt K ' 



d 2 w 
dz 2 



m 2 wo 



(8) 



d 2 w 

dr 2 
1 d 



3 dw 
r dr 



d 2 w 



+ m w 



+ --^-d ■ w(x) = 



dr 



(9) 



1 dw z 



dt 2 



r dr 



dr 2 
d 

—d ■ w(x) = 
dz 



dz 2 



+ m w z 



(10) 



using d ■ w(x) given in Eq. (J7JI. It is these equations 
which we solve to study the evolution of the IF-fields in 



the case of two colliding electroweak bubbles. For the 
calculations presented in this paper, the equations are 
solved in Mathematica [l2j using the built-in NDSolve 
function. 



Initial Conditions for Bubble Nucleation 

The auxiliary condition is maintained for all time if 
it (and its time derivative) is satisfied at an initial time 
to. This leads to initial conditions, or constraints, among 
the fields and their time derivatives at to [Hfl. In addi- 
tion to w z (to, x), we are also free to choose initial values 
for Wo(to, x), w(to,x) and dw z {to, x)/dt; for simplicity, 
we take the latter to vanish. The auxiliary condition, 
Eq. ([3]), then gives, in cylindrical coordinates, 



dw (t , x) dw z (t , x) _ w z (t , x) dm(t x) 



dt 



dz 



= 2 



m(to,x) dz 



(11) 



The initial condition for the time derivative of Wq (to, x) 
follows directly from Eq. (Ill)) and the boundary con- 
ditions BCI and BCII on w z (t 0l x). These boundary 
conditions are satisfied with w z (to,x) proportional to 
±m(to,x) when the colliding bubbles are well separated, 
which in turn may be satisfied if 



dw z (t ,x) _ w z (t ,x) dm(t x) 



dz 



m(to,x) dz 



(12) 



Combining Eqs. 
takes the form 



(|llll2p , the time derivative of w at to 



dw Q (to,x) _ ^dw z (t ,x) 



dt 



dz 



(13) 



The profiles of w z (to,x) for BCI and BCII are shown 
in Fig. [Tj It is seen that the bubbles collide shortly after 
nucleation with initial radii of r n — 20 in units of the 
inverse W mass. Corrections are required for the slight 
overlap of the bubbles and for fact that other bubbles 
give rise to an average scalar field, but these are small 
effects 10] and our calculated magnetic field should be 
reasonably accurate without them. 

The value of uu z (to, x) at t = to is fixed here by normal- 
izing the W fields inside the bubble as in Ref. Q to give 
a reasonable number of W^-bosons inside the bubbles 
under certain assumptions about the thermal conditions. 
By choosing the normalization in this way in both calcu- 
lations, a direct comparison with the results of Ref. § 
becomes more meaningful. 



Functional Form of the Scalar Field 



To solve the EOM in Eqs. (|8|9ll0|) . we need the func- 
tional form of mix) as a function of time. For a single 
bubble, m(x) is proportional to the magnitude of the 
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FIG. 1: Profile of initial conditions for w z 



scalar field p(x), which at zero temperature is simply the 
analytic continuation of the bounce solution to the Cole- 
man equation [l3| 



FIG. 2: Parameterized form for the scalar field of two colliding 
bubbles expending about their own centers and coalescing. 
The field is shown along the z-axis for the first few time steps 
for bubbles larger than the ones we have used to calculate the 
magnetic field. 



Current and Magnetic Field 



d 2 P (x) + p(x) 



dV 
dp{x) 



(14) 



in Euclidean space, where the bubble walls expand 
smoothly and retain the functional form of the bounce. 
At nucleation, we assume the scalar field is a simplified 
version of the bounce solution given in Eq. (6) of Ref. [| 
with 77VA = 2/3, 77 = 1, and A = 4/9; this corresponds 
to a bubble surface that falls from its 10% to 90% values 
over a distance of approximately 4.4. In this paper, the 
speed of the bubble walls is taken to be c, as might be 
expected for a very strong phase transition. It is known 
from previous studies of the electroweak phase transition 
that before the bubbles collide the walls reach a constant 
speed, where friction from the plasma and pressure in- 
side the bubbles balance, and that the bubble wall speed 
is definitely less than c [14]. It was found in that, 
for the case of Abelian bubble collisions, the resulting 
magnetic fields decrease in strength with decreasing wall 
speed, and we expect our fields to scale similarly. This 
issue was discussed further in 0] , and we plan to develop 
future calculations with a realistic wall speed. We as- 
sume that the bubbles nucleate at T = Tc = 166 GeV as 
in Ref. [|. 

For two colliding bubbles, rather than solving (|14p di- 
rectly, we choose a simple parameterized form for the 
scalar fields such that the bubble walls will expand 
smoothly, and that the scalar field will remain exactly 
constant throughout the bubbles in the broken symme- 
try phase. Graphs of this parameterization are shown in 
FigJU For calculating m(x) from the scalar field we take 
into account that there is also an average scalar field from 
the other bubbles in the medium, which in this work is 
assumed to have a magnitude of 10 % of the scalar field 
in the interior of a bubble. 



In cylindrical coordinates, using the form as before 

(15) 

the current |0J) can be written as [y = 0): 
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3x ^=3- 
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dt 
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r dr 
dw b 



r 
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2 <^r 



dr 



+w r- 



dr 



2wt 



dw b 



and for [y = 3): 



4ir Jz = ke abs [w b z (d ■ w a ) - 



dw ob 
dz 



+r 2 w a ^- + 2w a ^ 
oz at 

dw dw 

or oz 

Having obtained the current, we can now calculate the 
magnetic fields directly from the Maxwell equations 



dw b 



,dw b 



(16) 
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(18) 
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which may be written in terms of the vector potential 
as 



(20) 



(21) 



(22) 



(23) 



Taking the form of the vector potential in cylindrical 
coordinates to be 



d 2 A v -d u (d»A„) = j v 
Working in the axial gauge 

A z = 

the {y = z) equation becomes 

-d z {d^)=j z , 
which specifies the divergence of A^ as 

W = - / jzdz' . 



A, 



ao, f=0; 
ax v is=l,2; 
a z v=3 , 



(24) 



the equation for the vector potential a(x) becomes 



Numerical results and comparison to previous work 

Here we show the magnetic field with the boundary 
conditions in Fig. [1] and compare it to analogous results 
calculated from the theory of Ref . [8j , plotting the fields 
at comparable intervals 5t following the onset of the col- 
lision at to = 0. 
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FIG. 3: Magnetic field calculated in this work in the trans- 
verse direction at z=0 for times t =5 (solid curve), 10 (short 
dash curve), 15 (medium dash curve) and, 20 (long dash 
curve) in units where the w-boson mass mw = 1. The mag- 
netic field can be seen to be moving away from r = 0, and 
increasing in magnitude, as t increases. 
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r dr 



Jz dz')=j (25) 



Due to the axial symmetry of the collision, it can be 
easily shown that the only nonvanishing component of 
the magnetic field is B t f >1 given by 



da 

dz 



(26) 



The quantity may be found directly from Maxwell's 
equations by taking the derivative of (f2"5)l with respect to 
z and using equation ([26]) . 



cP__ d^l 
dt 2 dr 2 r 



4tt 



dr 
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dz 



dz 2 



(27) 



As a preliminary result, we solve (|2T|) on a coarse grid 
using an interpolating function for the current. Equa- 
tions (|8I9I10P were solved for t = to 30, r = .01 to 60 
and z — —80 to 80, and values for the current were cal- 
culated on a grid with a step size of 2 in the t, r and z 
directions. Mathematica [f| was then used to construct 
a polynomial interpolating function used as the current 
inEq. ((23). 



The corresponding fields calculated from Ref. [8j are 
shown in Fig. [U Bubble surface dynamics seems to pro- 
duce fields somewhat larger in magnitude, and it is ex- 
pected that bubble walls of even smaller surface thick- 
ness might grow even larger llj . The field in Fig. [3] can 
be seen to be more concentrated near the center of the 
bubbles, and for this reason will have a smaller scale at 
the completion of the phase transition. However, since 
the rate at which w z expands relative to the scalar field 
depends on the choice of dw z (to)/dt, it would be inter- 
esting to quantify in future work the extent to which the 
scale and magnitude of the magnetic field might increase 
with a choice different from the one made on this paper, 
dw z (to)/9t = 0. 

We have not attempted to determine the present day 
magnetic fields that are seeded by our fields generated 
during the EWPT since this is a complicated problem 
of plasma physics that has been studied extensively else- 
where. The most recent of these [lj| show the importance 
of helicity and supports the possibility that galactic clus- 
ter magnetic fields may be entirely primordial in origin. 



SUMMARY AND CONCLUSION 

We have shown how bubble surface dynamics affect 
magnetic seed field creation in collisions of bubbles in 
primordial first-order electroweak phase transition by ex- 
tending the study in Ref. § to treat, for the first time, 
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in physics because of its fundamental importance. 
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FIG. 4: Magnetic field calculated in [8J in the transverse di- 
rection at z—0 for times t —5 (solid curve), 10 (short dash 
curve), 15 (medium dash curve) and, 20 (long dash curve) in 
units where the w-boson mass mw = L 



the case of collisions of bubbles with walls of finite thick- 
ness. By working in the linear regime of gentle collisions, 
we are able to decouple the coupled and highly nonlin- 
ear partial differential equations that describe the evolu- 
tion of the scalar and w fields equations appearing in the 
Lagrangian of an appropriate extension of the Standard 
Model. This simplifies the equations, making the influ- 
ence of the surface dynamics relatively easy to study nu- 
merically. We find results in qualitative agreement with 
the previous studies, but allowing for the possibility that 
the magnetic seed fields could be even larger for suffi- 
ciently weak first-order phase transitions in which bub- 
bles occur with walls of even smaller thickness. Our work 
thus leaves open the possibility that the magnitude of the 
magnetic seed fields could be even larger than those that 
we find, making them an even more likely candidate for 
the origin of observed galactic and extra-galactic mag- 
netic fields. In this case the observation of these present 
day fields hold important clues to the form of the exten- 
sion of the Standard Model, a subject of intense interest 
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